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We present the analytical approach to scalar field theory on the fuzzy sphere 
which has been developed in hep-th/0706.2493. This approach is based on con- 
sidering a perturbative expansion of the kinetic term in the partition function. 
After truncating this expansion at second order, one arrives at a multitrace 
matrix model, which allows for an application of the saddle-point method. 
The results are in agreement with the numerical findings in the literature. 

x. Introduction 

It seems quite natural to expect that as one approaches the Planck scale, one has to re- 
place the smooth structure of spacetime by some form of quantized geometry. The usual 
quantization procedure as well as string theory suggest that the first modifications that 
should be encountered are noncommutative geometries with constant deformation param- 
eters. Besides the well-known Moyal-plane Rj and its 2(i-dimensional generalizations with 
their quantized function algebras based on = 18^, the so-called fuzzy geometries 

. have received more and more attention recently. 

A fuzzy geometry is essentially a noncommutative deformation of a Riemannian man- 
ifold which come with a Laplace operator with a discrete spectrum. Essentially, one 
truncates this spectrum and deforms the algebra of the corresponding truncated set of 
eigenfunctions to achieve closure under multiplication. The most prominent example of 
such a space is the fuzzy sphere [2] with a quantized function algebra based on the relation 



X 
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Since fuzzy spaces are described by function algebras with finitely many degrees of 
freedom, they might prove useful as regulators for quantum field theories. After calculat- 
ing the path integral depending explicitly on the deformation parameter 9, one should be 
able to recover the commutative path integral in a certain limit 9 — > 0. This regularization 
procedure would have several advantages over the lattice approach as, for example, it pre- 
serves a number of continuous symmetries and is not expected to suffer from the fermion 
doubling problem. Furthermore, numerical simulations are easily performed within this 
framework. 

A possible obstacle to using fuzzy geometries as regulators has been pointed out by 
various authors: the naive application of this procedure yields the wrong commutative 
limits in the case of fuzzy scalar field theory [3], as is also suggested by numerical studies 
of the phase diagram [I]. However, modifications of the action of the theory have been 
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proposed [5] and for scrutinizing them, an analytical handle on the partition function of 
this theory is desirable. In [I], such an approach was developed and we will review the 
approach and summarize the findings in the following. 

2. The fuzzy sphere and further fuzzy geometries 

To begin, let us briefly review the construction of a number of fuzzy geometries. First, 
recall that the complex space C" = R 2n is rendered noncommutative by replacing the 
complex coordinates z a , Zp with the creation and annihilation operators of n harmonic 
oscillators satisfying the algebra [a a , a^] = 5 a p. Functions become thus linear operators 
on the infinite-dimensional Fock space T generated by the creation operators from the 
vacuum |0). By normalizing the coordinates z a — > z a /\z\, we descend from C n to the 
sphere S 2n ~ x . To descend further to CP n ~ , we can use the Hopf fibration 

1 -> U(l) -> S 2 "" 1 -> CP"" 1 -> 1 , (2.1) 

which tells us that functions on CP n_1 are obtained from the functions on g 2 ™ -1 by 
factoring out a U(l)-action. We can take this action to be z a \— > e 1(p z a , and thus we see 
that functions on CP™ -1 are built from monomials containing an equal number of z a and 
z a . As a basis for these functions, we use 

z ai ... z aL z f3l ...zp L , (2.2) 

which under the above quantization prescription turns into the operator basis [B] 

al i ...ai L \0)(0\ap 1 ...a PL . (2.3) 

This basis spans the space of linear operators acting in the L-particle Hilbert space Hi 
of the above Fock space T. We have dim(7i £ ) = ("~ 1 )-( Tt+z ') ) ari d i n particular, for the 
sphere S 2 = CP 1 , dim(7-^) = L+ 1. Functions on the fuzzy sphere can thus be encoded in 
(L+ l) 2 -dimensional matrices; real functions correspond to hermitian matrices. Note that 
this map between monomials and operators also gives a direct quantization prescription 
for the function algebra of a projective algebraic variety as discussed in [TJ. 

Note also that the quantization of CP n as described above corresponds to the geo- 
metric quantization (Toeplitz quantization) of CP™ with the quantum line bundle Jzf = 
0(1)® = O(L): Recall that this quantization procedure consists of replacing the alge- 
bra of smooth functions on a projective algebraic variety, C°°(M), with the endomor- 
phisms of sections End (T(M, Jzf)) of the quantum line bundle =5f. As the set of sections 
T(CP 1 , 0(L)) is spanned by z ai . . . z aL , the quantized algebra is spanned by the operators 
Zai ■ ■ ■ z a L ~d^~ ■ ■ ■ d^T an d ^ i s thus equivalent to the operator basis (12.31) . 

To capture the topology of the space we are quantizing, it is crucial to specify an 
additional structure, which is usually taken to be a Dirac operator in noncommutative 
geometry. For our purposes, however, it is sufficient to define a Laplace operator. Its 
definition is most easily gleaned from an alternative, group theoretic point of view. 

Consider a finite dimensional, irreducible representation p of SU(ra), extended to a rep- 
resentation of U(n). Such a representation can be labeled by a Dynkin diagram endowed 
with Dynkin labels d, . . . , a n -i'- 

o o o 

U(l) U(l) U(l) U(l) U(l) (2.4) 
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where we have arranged the Cartan generators of the maximal torus U(l) xn around the 
corresponding roots <5j corresponding to the Dynkin labels aj. To every simple root, we 
have a pair of raising and lowering operators, E s ., and the Dynkin label indicates the 
highest non-trivial action of the lowering operator on the highest weight state \p) of this 
representation: 

(E^riri * o = (£ Q - s r+v> • (2.5) 

Putting Dynkin labels to zero enlarges the isotropy group of the highest weight state 
in p, as the generators E s , now leave \p) invariant and combine the neighboring 
U(n) x U(m) in the diagram (12 4[) to U(n + m). One can now show that there is a one- 
to-one correspondence between points p on the coset space U(ri)/(U(m 1 ) x . . . x U(m fc )) 
and coherent states \p) in a representation p in which the isotropy group of the highest 
weight state \p) is U(mi) x . . . x U(mfc). 

An obvious quantization prescription is now to associate an operator / to a function 
f(p) according to 

m = (pirn , (2.6) 

where \p) is the coherent state in p corresponding to the point p. In this way, one can 
find quantizations for all complex flagmanifolds and their supersymmetric extensions as 
shown in [H]. 

Recall also that the Young diagram corresponding to the representation p(a\, . . . , a„_i) 
is given by 



n — 1- 



(2-7) 



which in the case of S 2 = CP 1 = U(2)/(U(1) x U(l)) just becomes p{a 1 ) with a x = L or 

= span(at i ...atjO) ) (2.8) 



A function is therefore mapped via (12.61) to a linear operator, 




- span(dt i ...atjO)(0|a /3l ...a /3L ) , (2.9) 

and we arrive again at the above quantization procedure. 

The representation of the generators of the isometries on CP 1 = SU(2)/U(1) on states 
in p is the usual Schwinger representation, = a^cLgO^. The Laplace operator is then 
the square of this action on functions, which is also the second Casimir operator in p: 

Af - A/ := C 2 f = [Li, [Li, /]] . (2.10) 

The requirement that the measure involved in the integration procedure is invariant under 
the isometries of CP 1 as well as the usual rule of partial integration (f M da = J dM a = 0, 
as the spaces M under consideration are compact) and the volume formula on the sphere 
force us to define integration according to f g2 dA f — > tr (/). 

We have now a complete quantization of the space CP 1 at hand, and we can easily 
translate a commutative scalar field theory to the noncommutative setting. From now on, 
we will omit hats over operators for convenience. Also, in accordance with the standard 
nomenclature of matrix models, we will label the size of the matrices encoding functions 
by N = L + l. 
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3. Fuzzy scalar field theory 



With the prescription of the last section, we immediately arrive at the following action 
for scalar 4 -theory on the fuzzy sphere: 

S = 7 tr (-^$C 2 $ + r$ 2 + s$ 4 ) = 7 tr $] [L,, $] + r $ 2 + ^ $ 4 ) , (3.1) 

where we introduced the shorthand notation 7 = For simplicity, we will put R = 1 

in the following. The partition function of the model reads as 

Z = J d MD ($) e"^ = J dfl D ($) e -^tr(-f[L l ,$][L^] + ,^ +9 ^) ^ (3 2) 

where the measure d/i£)($) denotes the Dyson measure on the set of hermitian matrices 
of dimension N x N, d/i D ($) := Ui<j d^(%) Ui<j dS^-). 

It is well-known that the pure matrix model limit a = of this model has a third 
order phase transition in the large N limit at g = j^r 2 = -§^t 2 for r < 0. At this 
parabola in the left half of the r-g-plane, the double well potential becomes sufficiently 
deep for the support of the eigenvalue density to split into two disjoint pieces. 

Numerical studies on the lattice [9] of the planar commutative limit of this model, 

^e"/ d SW)W ; (3.3) 

indicate a second-order phase transition for g/r ~ —10.24, i.e. at a line in the r-g-plane. 

Interestingly, numerical studies [I] of the fuzzy scalar field theory (13. ip find a com- 
bination of both phase transitions with a triple point at which the parabola meets the 
line. This seems at least to suggest a contradiction to the idea that the fuzzy sphere can 
be used as a regulator for scalar quantum field theory on the plane. After gaining an 
analytical handle on the phase diagram, we will therefore apply our technique to study 
modifications proposed to cure this problem and we will see that the triple point is shifted 
off to infinity; this is an indication for an improved situation as the phase diagram turns 
into the one of (13. 3p . 

Let us briefly recall why scalar field theory on the fuzzy sphere is more difficult than 
the more common matrix models. First, let us consider again the matrix model limit 
a = 0, which is exactly solvable [TU]. Using the decomposition $ = fiAfit, where A = 
diag(Ai, . . . , Ajv) and Q is a unitary matrix, the dependence on Q drops out from the 
action due to cyclicity of the trace. The partition function (I3.2p turns into the one of the 
eigenvalue model 



Z a=0 = J $)\ A 2 (A) J dfi H {Q) e-'W'-^+'EA*) 

= J ^ e - 2 E,> 3 ta|A.-A J |-/37(rE,^+sE 1 Af) ; 



(3.4) 



where 3>\ := ]^[ i=1 dAj, d/iff(f2) is the Haar measure on U(iV) and A (A) is the Van- 
dermonde determinant A(A) := det([A^~ 1 ]jj) = rij>j(^« — From here, one can 
continue with various techniques to evaluate Z, as, e.g., the saddle-point approximation 
or the method of orthogonal polynomials. The latter even yields an exact result for Z at 
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finite N. However, f2 does not commute with the external matrices Lj appearing in our 
model (13.11) and thus this method is not directly applicable in our case. 

More generally, there is a solution for hermitian matrix models of the form S = 
tr (Vi(A<&) + V2($)) with a single external matrix A as shown in [TT]: One can use a 
character expansion together with the orthogonality relation for elements of U(iV) in an 
irreducible representations p, 

J dp H (n) Xp (Art^) = j-±^ Xp (A) Xp (A) , (3-5) 

to arrive at a closed formula for the partition function. Unfortunately, having three 
external matrices as in our action (13.11) renders the character expansion so complicated 
that it is essentially useless for our purposes. 



4. Calculating the phase diagram 

Before performing the perturbative calculation, we can employ symmetry arguments 
to make some statements about the expected results. First of all, we recall that the 
Dyson measure is invariant under the adjoint action by a unitary matrix Q: dpoi^) = 
dpolQ^Q^). This implies that we can replace the action (13.11) . which is not invariant 
under this action, by an effective action S e g with this invariance under the functional 
integral: 

J dp D (®)e- s = J d/i D ($)e-^ , e- s -[*] = 1 J dp H (n) e~ s ^ . (4.1) 
Because of this invariance, the effective action has to be of the form 

s cS = s « tr (* n ) + E tr tr ( $m ) + E s ™ k tr tr ( $m ) tr ( $fc ) + • • • ' 

n n,m n,m,k 

(4.2) 

which can now be trivially recast into an eigenvalue model S , e fr(A). 

Furthermore, symmetry arguments severely restrict the form in which the kinetic term 
can appear in the effective action S e g(A). Since C 2 t = 0, the kinetic term tr ($(72$) 
can only depend on the difference of eigenvalues Aj — Xj. Because of the Z^-symmetry 
Aj — > —A», we have to take this difference to an even power and the permutation symmetry 
between the eigenvalues demands that we sum over all i > j. Altogether we arrive at 

Seff(A) - V(A) = ^(mi,ni)...(m fc ,n fc ) S 2mi • • • S 2m fc > ( 4 -3) 

fe,mi,ni,...,mfe,nfc 



where 

^2m := (2 2m ) n := (E( A ^ A i) 2m ) " ( 4 ' 4 ) 

i>j 

The idea of treating the kinetic term perturbatively is known as the hopping parameter 
or high-temperature expansion and this technique has been successfully used to analyze 
two-dimensional scalar 4 -theory on the lattice, see e.g. [12] and references therein. In our 
case, it will allow us to treat the partition function in principle exactly for any truncation 
of the perturbative series. 
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We will restrict our considerations to an expansion of the kinetic term up to second 
order: 

o2 2 2 

e ^a^K ab ^ = x + p ia $* Kab <$>h + tLlJL<$><>K ab $ b § c K cd § d + C(a 3 ) , (4.5) 

2 

where 

K ab = tr{[L h r a ][L h r b \) and $ = $% = + $<V a (4.6) 

and t q are the Gell-Mann matrices of SU(iV) normalized according to tr (r a r b ) = 8 ab . To 
integrate over the angular variables, we have to compute 

J dii H {fyK ab tr{T a nArt)tr{T b nArt) with $ a = tr^fiA^). (4.7) 

Having in mind the orthogonality relation ( 13. 5ft . we rewrite 

tr ((r^A^) ® (r^Afit)) = tr ((r a ® r 6 )(fi ® fl)(A ® A)(fi f ® fi f )) . (4.8) 

Splitting the tensor product □ <g> □ into irreducible representations p, we arrive at the 
formula 

/ dfi H (n)K ab ^ b = K ab ^ dht(p) tr p( r<1 ® r& ) tr p(A <8> A) , (4.9) 
^ p 
where tr p denotes the trace in the irreducible representation specified by p. The explicit 
expressions for tr p in all relevant representations are found in pQ, and using these, one 
straightforwardly computes that, in accordance with ( 14.31) . 

f dv H ($)K*&* b = f E( A *- A i) 2 = / d ^( $ )^ 2 ( tr ( $2 )-^ tr ( $ ) 2 ) • 

At second order, we obtain after a more tedious but also straightforward calculation 
a more complicated result, which is again of the form predicted in (14.31) . Instead of pre- 
senting it here, let us directly jump to the full large N limit of the model. For this, we 
re-exponentiate the contributions found perturbatively into the effective action. Neglect- 
ing terms which are subdominant in N, we arrive at the following eigenvalue model: 



S = 7 ^(rA, 2 + ^)+ 7 ^(-fiV(A J -A J ) 2 + ^iV 2 (A J -A J ) 4 -fln|A l 



A, 



The transition to continuous variables is performed as usual by rescaling Aj — > A(^) = 
X(x) with < x < 1 and turning the sums into integrals: J^Li — > A" dx. Moreover, a 
common power of A" has to be factored out from every term in the action. This power is 
determined by the logarithmic term to be A^ 2 and yields a rescaling, 

a = N 9a a , r = N 9r f , g = N e <>g and X(x) = N 9x X(x) , (4.10) 

with values for 9 g and 9 r which remarkably are found to be consistent with the 
ones obtained numerically in [4]. Altogether, we have now the partition function Z = 
J @X exp(-N 2 S) with 

1 r 1 



S = An dx(fX 2 (x)+gX\x) + / dyl - j(X(x) — X(y)Y 

Jo V (4.11) 

+ *f{\(x) - X(y)r - i In |A» - A(y)|)) , 



which we can evaluate using the saddle point approximation. For this, we introduce 
as usual the eigenvalue density u(A) = dx/dX and follow the canonical procedure for 
determining tt(A) (see e.g. [ID] , cf. [13]), which yields in the single cut regime 

u(A) = Uf-a + I2ix~a 2 c 2 + 4 (g + 5 2 + 8 (g + *f ) A 2 ) V5 2 - A 2 . (4.12) 

Here, C2, the second moment of w(A), is determined by a self-consistency condition. This 
solution enables us to locate a phase transition, i.e. a curve C in the r-g-plane, at which 
the eigenvalue density u(X) becomes negative. This curve is found to be given by 

= ^ (-63S 2 + 16f 2 ± (4f - a)Vl6f 2 - 8af - 955 2 - 8ar) . (4.13) 

Various arguments suggest to identify the turning point of this curve with the triple point 
of the phase diagram predicted numerically. The location of this point for a — 1, (r, g) « 
(—2.7, 0.25), is in good agreement with the numerical result (r, g) = (— 2.3±0.2, 0.52±0.02) 
and justifies further this identification. 

As mentioned in the introduction, modifications to the action of fuzzy scalar field 
theory are necessary, if one wants to employ fuzzy spaces to regularize 4 -theory on the 
plane. The modification proposed in [5j takes into account a wave function renormalization 
and - in its simplest form - reads as 

S = 7 tr (a$(C 2 + /sAC , 2 )$ + r$ 2 + £$ 4 ) . (4.14) 

This modification is readily treated in our formalism, as it amounts to 

K ab -> K ab := K ab + KK ac K cb and a -> a = a(l + §«) . (4.15) 

Increasing a now moves the turning point of C - and thus the triple point - off to infinity. 
We thus confirmed that introducing k has the desired effect. 



5. Summary and future directions 

Altogether, we achieved the following: We formulated a generalized character expansion 
technique for the treatment of fuzzy scalar field theory. This technique yields exact 
results at any order in a perturbative expansion of the kinetic term and can be evaluated 
- in principle - straightforwardly. It is worth stressing that this approach is directly 
applicable to field theories on the fuzzy sphere with more general potential and it is 
readily adapted to scalar field theories on other fuzzy spaces. We used this expansion 
technique to reformulate scalar field theory on the fuzzy sphere as a multitrace matrix 
model. The results of the approximation look promising and motivate further studies. 

An important task for future work will be to explore the full set of one- and two-cut 
solutions as well as their role in the explanation of the phase diagram. The qualitative 
effects of higher order corrections have to be taken into account, if the full analysis at 
second order should not suffice to explain the phase structure. Obviously, our technique 
should also be applied to the probably more interesting case of four dimensional theories. 
Finally, one might wonder whether there is a connection between the multitrace matrix 
model we found from fuzzy scalar field theory and the original motivation for studying 
multitrace matrix models [13j : the definition of c > 1 string theories. 
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